Fluctuation Probes of Quark Deconfinement 



Masayuki Asakawa 1 , Ulrich Heinz 2 , and Bcrndt Miiller 3 
1 Department of Physics, Nagoya University, Nagoya 464-8602, Japan 
2 Theoretical Physics Division, CERN, CH-1211 Geneva 23, Switzerland 
3 Department of Physics, Duke University, Durham, NC 27708-0305 
(final version, February 1, 2008) 



o 
o 
o 

(N 



ON 

m 
> 

on ; 
\o ■ 

m : 
o . 
o ■ 
o ; 

^ : 

i 

Q . 
^ ■ 



X 



The size of the average fluctuations of net baryon num- 
ber and electric charge in a finite volume of hadronic matter 
differs widely between the confined and deconfined phases. 
These differences may be exploited as indicators of the for- 
mation of a quark-gluon plasma in relativistic heavy-ion col- 
lisions, because fluctuations created in the initial state survive 
until freeze-out due to the rapid expansion of the hot fireball. 



Fluctuations in the multiplicities and momentum dis- 
tributions of particles emitted in relativistic heavy-ion 
collisions have been widely considered as probes of ther- 
malization and the statistical nature of particle produc- 
tion in such reactions Jl| The characteristic behav- 
ior of temperature and pion multiplicity fluctuations in 
the final state has been proposed as a tool for the mea- 
surement of the specific heat and, specifically, for the 
detection of a critical point in the nuclear matter phase 
diagram || . Although the hot and dense matter created 
in heavy-ion collisions is not directly observed at the cri- 
tical point (if one exists) but rather at the point of ther- 
mal freeze-out where particles decouple from the system, 
certain features of the critical fluctuations were shown to 
survive due to the finite cooling rate of the fireball Q . 

We here draw attention to a different type of fluctu- 
ations which are sensitive to the microscopic structure 
of the dense matter. If the expansion is too fast for lo- 
cal fluctuations to follow the mean thermodynamic evo- 
lution of the system, it makes sense to consider fluctu- 
ations of locally conserved quantities that show a dis- 
tinctly different behavior in a hadron gas (HG) and a 
quark-gluon plasma (QGP). Characteristic features of 
the plasma phase may then survive in the finally observed 
fluctuations. This is most likely if subvolumes are consid- 
ered which recede rapidly from each other due to a strong 
differential collective flow pattern as it is known to exist 
in the final stages of a relativistic heavy-ion reaction. 

Three observables satisfy these constraints and are, 
in principle, measurable: the net baryon number, the 
net electric charge, and the net strangeness. Here we 
will focus on the first two as probes of the transition 
from hadronic matter to a deconfined QGP. Because they 
are sensitive to the microscopic structure of the matter, 
their unusual behavior would provide specific informa- 
tion about the structural change occurring as quarks are 



liberated and chiral symmetry is restored at high tem- 
perature. Our proposal differs from recent suggestions 
involving fluctuations in the abundance ratios of charged 
particles [ [lO[ and in the baryon number multiplicity pd| ] 
in that we only consider locally conserved quantities. We 
also disregard dynamical fluctuations of the baryon den- 
sity caused by supercooling and bubble formation fll^] . 

We consider matter which is meson-dominated, i.e. 
whose baryonic chemical potential /i and temperature T 
satisfy \x < T. Our arguments will thus apply to heavy- 
ion collisions at CERN SPS energies and above. In the 
following, we first explain qualitatively how hadronic and 
quark matter differ with respect to net baryon number 
and electric charge fluctuations. We then present ana- 
lytical calculations supporting the argument. Finally, we 
estimate the rate at which initial state fluctuations are 
washed out during the expansion of the hot matter in the 
final, hadronic stage before thermal freeze-out. 

In a hadron gas nearly two thirds of the hadrons (for 
/i <C T mostly pions) carry electric charge ±1. In the de- 
confined QGP phase the charged quarks and antiquarks 
make up only about half the degrees of freedom, with 
charges of only ±| or ±|. Consequently, the fluctua- 
tion of one charged particle in or out of the considered 
subvolume produces a larger mean square fluctuation of 
the net electric charge if the system is in the HG phase. 
For baryon number fluctuations the situation is less ob- 
vious because in the HG baryon charge is now only car- 
ried by the heavy and less abundant baryons and an- 
tibaryons. Still, all of them carry unit baryon charge ±1 
while the quarks and antiquarks in the QGP only have 
baryon number ±i. It turns out that, as /i/T — ► 0, 
the fluctuations are again larger in the HG, albeit by a 
smaller margin than for charge fluctuations. At SPS en- 
ergies and below the difference between the two phases 
increases since the stopped net baryons from the incom- 
ing nuclei contribute to the fluctuations, and more so in 
the HG than in the QGP phase. 

Generally, if O is conserved and \i is the associated 
chemical potential, in thermal equilibrium the mean 
square deviation of O is given by 
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where (O) = Tr O e -CH-»o)/T ^ ( 

Nb the r.h.s. of (Q) is T times the baryon number sus- 
ceptibility which was discussed earlier in the context of 
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possible signatures for chiral symmetry restoration in the 
hadron-quark transition [T^ ] . 

In general, the relative fluctuation of any extensive 
variable vanishes in the thermodynamic limit V — > oo be- 
cause the expectation value (O) increases linearly with 
the volume V while the fluctuation AO grows only like 
VY. In reality, the value of a conserved quantum num- 
ber of an isolated system does not fluctuate at all. How- 
ever, if we consider a small part of the system, which is 
large enough to neglect quantum fluctuations, but small 
enough that the entire system can be treated as a heat 
bath, Eq. (g) can be used to calculate the statistical un- 
certainty of the value of the observable in the subsystem. 
This is the scenario considered here. 

We first discuss the fluctuations of the net baryon num- 
ber. Since baryons are heavy, in the dilute HG phase we 
can apply the Boltzmann approximation |L4j : 



N±(T, fj.) = A±(T, 0) exp(± M /T) . 



Here Njj 1 denotes the number of baryons (+) and an- 
tibaryons (— ), respectively. The net baryon number is 
Then the net baryon number fluctua- 
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tions in the hadronic gas are given by 

(AN b ) 2 HG = N+ + N b -=2 N±(T, 0) cosh( M /T) . (3) 

This result makes sense, because the fluctuation of cither 
a baryon or an antibaryon into or out of the subvolume 
changes the net baryon number contained in it. 

To estimate (AN b ) 2 in the QGP phase, we use the 
exact result for the baryon number density in an ideal gas 
of massless quarks and gluons (for two massless flavors): 
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where V denotes the volume of the considered subsystem. 
It is convenient to normalize this by the entropy density 
(again for two quark flavors plus gluons): 
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The later expansion being nearly isentropic, the ratio 
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provides a useful measure for the fluctuations predicted 
for a transient quark phase. The entropy can be estima- 
ted from the final hadron multiplicity Jl6|] . 

For high collision energies (fi/T— >0), the ratio (@) ap- 
proaches a constant; even for SPS energies, the //-de- 
pendent correction is at most 5%. The many resonance 
contributions make it difficult to write down an analytic 
expression like (j^) for the entropy density in a hadron 
gas, but it is clear that the stronger //-dependence of (0) 



compared to (g) induces a stronger //-dependence of the 
corresponding ratio @ in the HG phase. This translates 
into a stronger beam energy dependence of the ratio (|^) 
near midrapidity in the HG than in the QGP phase. 

Before providing numerical illustrations, let us com- 
pare these results with those for net charge fluctuations. 
All stable charged hadrons have unit electric charge; 
again using the Boltzmann approximation, which only 
for pions introduces a small error of at most 10%, we 
find 
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where 7V c h is the total number of charged particles emit- 
ted from the subvolume. To find the expression for a 
noninteracting QGP, we introduce the electrochemical 
potential cj> which couples to the electric charges q u = | 
and qd = — | of the up- and down-quarks: 



(2) (g_W) _ 
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We differentiate with respect to <f) at <j> = and normalize 
to the entropy density: 
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This is a factor | larger than the corresponding ratio (|^) 
for baryon number fluctuations, due to the larger elec- 
tric charge of the up-quarks, but shows the same weak 
//-dependence. The main difference to baryon number 
fluctuations arises in the HG phase: Since at SPS and 
higher energies the r.h.s. of (g) is dominated by pions and 
meson resonances, its //-dependence is now also weak. In 
contrast to baryon number fluctuations, charge fluctua- 
tions thus show a weak beam energy dependence in either 
phase, and only their absolute values differ JlTj . 

We now give some numerical values for the fluctu- 
ation/entropy ratios at SPS and RHIC/LHC. At the 
SPS, the net baryon number per unit of rapidity is mea- 
sured: d Nb /dy ~ 92 Jl8[ . The antibaryon/baryon ratio is 
« 0.085 JlJJlij], corresponding to dN b ~ /dy « 8.5. Com- 
bined with a specific entropy of S/N b 36 p0| , Eq. (||) 
thus gives (AN b ) 2 / S ~ 0.033 if the fluctuations reflect an 
equilibrium HG. If they have a QGP origin, Eq. (|^) gives 
(AN b ) 2 /S w 0.014 |J, i.e. about a factor 2.4 less. - 
The charge fluctuations in a HG can be evaluated from 
the measured charged multiplicity density at midrapid- 
ity, dN ch /dy « 400 @^9@, after correcting for res- 
onance decays [ fL6| . Assuming hadrochemical freeze-out 
at T w 170 MeV & 60% of the observed pions stem 
from such decays [§2|. One finds (AQ) 2 /5 w 0.06. If 
the charge fluctuations arise from a QGP, Eq. (§) gives 
{AQ) 2 /S re 0.036, i.e. 60% of the HG value. 

It is instructive to extrapolate these results to RHIC/ 
LHC energies (i.e. n/T— >0). We again assume hadro- 
chemical freeze-out at 170 MeV and use the particle 
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multiplicities predicted by hadrochemical models |^3|,^| . 
One obtains (AN b ) 2 /S « 0.020 in the HG, compared to 
0.0137 in the QGP, and {AQf / S « 0.067 in the HG 
phase, compared to 0.034 in the QGP. Only the first 
of these four numbers, corresponding to the hadronic 
baryon number fluctuations, changes by more than 10% 
as one proceeds from SPS to RHIC (see Fig. 0). 

These estimates, including our corrections for reso- 
nance decays, refer to ideal gases in equilibrium. Future 
work should address interaction effects on the thermal 
fluctuations in HG and QGP and treat resonance decays 
kinetically. We also point out potentially important non- 
equilibrium aspects: The fluctuation/entropy ratios in 
the QGP will be even lower (facilitating the discrimina- 
tion against HG) if initially the QGP is strongly gluon- 
dominated [p5| and hadronizes before the concentrations 
of the (baryon) charge carriers q,q saturate |2(|, or if 
hadronization itself generates additional entropy. 

We now discuss whether the difference between the two 
phases (typically a factor 2) is really observable. Even 
if a QGP is temporarily created in a heavy-ion collision, 
all hadrons are emitted after re-hadronization. Thus, it 
is natural to ask whether the fluctuations will not always 
reflect the hadronic nature of the emitting environment. 
We must show that the time scale for the dissipation of an 
initial state fluctuation is larger than the duration from 
hadronization to final particle freeze-out. It is essential 
to our argument that fluctuations of conserved quantum 
numbers can only be changed by particle transport and 
thus are likely to be frozen in at an early stage, similar 
to the abundances of strange hadrons, which are frozen 
early in the reaction and may even reflect the chemical 
composition of a deconfined plasma |27| . 

For our estimate we assume for simplicity that the fire- 
ball expands mostly longitudinally, with a boost-invarant 
(Bjorken) flow profile. Longitudinal position and rapid- 
ity are then directly related. Strong longitudinal flow 
exists in collisions at the SPS [p8[ , and the Bjorken pic- 
ture is widely expected to hold for collisions at RHIC 
and LHC. Consider a slice of matter spanning a rapidity 
interval Arj at the initial time Tj. (r is the proper time 
and r/= tanh _1 (z/i).) Its proper volume is Vj — AriArj 
where A is the transverse area of the fireball. We denote 
the initial total baryon density by pi = Pb+j>( T i)- We as- 
sume Ti = 170 MeV and T f = 120 Me V for the initial and 
final temperature p9| , corresponding to Tj w 2.5 fm/c and 
Tf w 7fm/c at the SPS, and n ~ 5 fm/c and r/ w 14 fm/c 
at RHIC. 

Let us first give a qualitative argument j3Q] for the 
survival of a baryon number fluctuation within a rapidity 
interval Arjfa 1. Between Tj and tj, this interval expands 
from a length of 5 fm to 14 fm (we use the RHIC num- 
bers here). Baryons have average thermal longitudinal 
velocity component v z =^v where v = (\v\) = ^&T/itM 
is the mean thermal velocity (v = 0.65 for baryons with 
M = 1 GeV at T= 170 MeV). Without rescattering, be- 



tween Ti and Tf a baryon which is initially at the center 
of this interval can travel on average only about 3 fm in 
the beam direction; hence it will not reach the edge of 
the interval before freeze-out. Because of rescattering in 
the hot hadronic matter, the baryon number actually dif- 
fuses more slowly, and a fluctuation will even survive in 
a smaller rapidity interval. 

For a quantitative argument, we need to estimate the 
flux of baryons in and out of the considered rapidity in- 
terval. Two effects need to be evaluated in this context. 
First, the difference in the baryon densities inside and 
outside the subvolume causes a difference in the values 
of the mean flux of baryons into and out of the volume. 
Denoting by v(t) the average thermal velocity of baryons, 
one finds that the initial fluctuation decays exponentially: 



AiV b (T) = AJV« exp 
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In the Bjorken scenario, the temperature T falls as t 1 ' 3 
and one finds for the remaining fluctuation at freeze-out 



AN b {T f ) = AJV« exp 



(Tf/T,) 1 / 2 ]) . (11) 



For the numbers considered here, the exponent is very 
close to —Vi/(2Arf), implying that the fluctuation sur- 
vives if Arj is larger than u,/2 w 0.33. 

The second effect that can wash out the initial fluctu- 
ation is fluctuations in the baryon fluxes exchanged with 
the neighboring subvolumes. These could eventually re- 
place the initial fluctuation with a thermal fluctuation 
that is characteristic of the conditions at freeze-out. The 
total number of baryons entering N^ en ' or leaving 
the subvolume between Tj and Tf is given by 

N (en) = N W = A £ Mr) - (r)dr _ (12) 

A similar calculation yields N^ n] = N^ lv) « N^v t /2A7j. 

N^ en ' and iV^ 1 ^ fluctuate independently; one therefore 
expects that the ratio of the mean square fluctuation of 



the number of exchanged baryons N b 
initial fluctuation is: 
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which is smaller than unity for A77 > x>i sa 0.65. 

We conclude that the short time between hadroniza- 
tion and final freeze-out precludes the readjustment of 
net baryon number fluctuations in rapidity bins Ar\ > 1 . 
A similar calculation applies to net charge fluctuations. 
Several refinements of our estimate are possible but are 
expected to partially cancel each other: Additional trans- 
verse expansion lets the temperature drop faster than in 
the Bjorken scenario. During hadronization cooling is 



3 



imp eded by the large change in the entropy density be- [11] 



twe qn QGP and HG. And finally, the short mean free 



path of baryons in hot hadronic matter will significantly 
reduce our above estimates of the dissipation of an initial 
state fluctuation. 

In conclusion, we have argued that the difference in 
magnitude of local fluctuations of the net baryon number 
and net electric charge between confined and deconfincd 
hadronic matter is partially frozen at an early stage in 
relativistic heavy-ion collisions. These fluctuations may 
thus be useful probes of the temporary formation of a de- 
confined state in such collisions. The event-by-event fluc- 
tuations of the two suggested observables for collisions 
with a fixed value of the transverse energy dEf /dy or of 
the energy measured in a zero-degree calorimeter would 
be appropriate observables that could test our predic- 
tions. Further discrimination can be achieved by measur- 
ing the beam energy dependence of the fluctuations: In 
the QGP the ratio (AQ/AN b ) 2 = | of charge to baryon 
number fluctuations is a beam-energy independent con- 
stant; in the HG phase it shows a significant beam energy 
dependence between SPS and RHIC/LHC energies. 
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Note added: After finishing this work we received a 
paper by Jeon and Koch plj who discuss similar issues. 
At the SPS they get (AQ) 2 /S'| HG w 0.13 which is more 
than twice our value, due to a smaller resonance decay 
correction to (AQ) 2 (30% instead of our 50%) and their 
omission of a 35% extra contribution to S |l6[ | from heavy 
particles (mostly the net baryons and strange hadrons). 
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FIG. 1. Schematic drawing of the beam energy depen- 
dence of the net baryon number and charge fluctuations per 
unit entropy for a hadronic gas and a quark-gluon plasma. 
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